Optics Letters

Minimally invasive surgery is a powerful technique that enables operations deep within the body while minimizing patient trauma and recovery time. Optical endoscopes are key to providing intraoperative vision but still face challenges due to the loss of essential senses, including depth perception and tactile feedback for tissue evaluation. Thus, it is critical to develop endoscopic imaging technologies that can augment operators with critical information. In this work, we explore a prototype multimodal 3D imaging endoscope that integrates volumetric light-field imaging with laser-speckle contrast imaging to simultaneously capture 3D structure and blood-flow information in a clinically relevant form factor.

Inverse-scattering methods enable label-free, quantitative visualization of a samples three-dimensional (3D) refractive index (RI), providing intrinsic and volumetric morphological contrast without exogenous labels. This is achieved by developing computational frameworks that reconstruct the samples 3D RI from a series of scattering measurements acquired under different data-capture conditions. Recent advances have demonstrated successful 3D RI reconstructions in multiple-scattering samples using angle-varying illuminations; however, these studies have primarily focused on non-absorptive samples. Here, we extend the multi-slice beam propagation (MSBP) inverse-scattering framework to reconstruct complex-valued RI, encompassing both the samples conventional RI (real part) and absorptivity (imaginary part). We show that reconstructing complex-valued RI makes the inverse problem ill-posed under angle-varying illumination alone, and that incorporating measurement diversity from both angle-varying illumination and sample defocus is necessary to ensure stable and accurate convergence. Experimental demonstrations were conducted on 1) dyed microsphere samples to characterize accuracy of reconstructed RI and absorptivity; and 2) diverse absorptive scattering samples to demonstrate biological utility. These results represent an important step for label-free volumetric imaging in biological tissue, which typically exhibits both scattering and absorption.

1.We report a repetition-controllable gain-managed nonlinear fiber amplifier (GMNA) that delivers near-infrared 50-fs pulses with pulse energies up to 150 nJ and a widely tunable repetition rate from 1-20 MHz, while maintaining stable pulse quality across the full range. Using this source, we demonstrate label-free multiphoton imaging--including metabolic autofluorescence (2PF/3PF), second/third-harmonic generation, and Simultaneous Label-free Autofluorescence Multiharmonic (SLAM) microscopy imaging--across live cells, human lung spheroids, and hard tissues. We further assess the impact of laser repetition rate on photodamage at fixed pulse energy, supported by preliminary measurements indicating lower damage at lower repetition rate. Collectively, the compact architecture and repetition-rate agility of the GMNA enable real-time optimization of imaging speed, depth, and sample safety for advanced biological microscopy.

This study presents the use of sub-micron nanobubbles (NBs) as contrast agents for ultrasound localization microscopy (ULM), a super-resolution imaging technique that visualizes microvascular structure and flow beyond the acoustic diffraction limit. While ULM has traditionally relied on micron-sized microbubbles (MBs), the reduced dimensions and prolonged circulation times of NBs make them attractive candidates for localization-based imaging. However, their weaker acoustic responses present significant challenges for reliable detection and tracking. To address this challenge, we developed the ULM Master GUI, an interactive framework for optimization of the complete ULM processing pipeline. Using custom ultrasound-compatible wall-less gelatin flow phantoms containing vessel-mimicking channels and bifurcations ranging from 100 to 500 m, we demonstrate that NB-based ULM achieves velocity reconstruction and flow partitioning measurements comparable to conventional MB-based ULM. Across all investigated geometries, NBs faithfully reproduced the underlying flow patterns and hemodynamic behavior despite their substantially reduced acoustic scattering. These findings establish the feasibility of NB-based ULM, expand the range of contrast agents available for localization microscopy, and provide a foundation for future super-resolution ultrasound imaging using nanoscale acoustic contrast agents. The ULM processing GUI is publicly available at https://github.com/grisha1998/ulm-super-resolution-toolbox.

AO_SCPLOWBSTRACTC_SCPLOWEarly tumor detection is critical for improving patient survival and recovery. Clinically, tissue palpation is routinely used to identify regions of abnormal stiffness, a hallmark of many pathological conditions. However, palpation is restricted to anatomically accessible sites and remains highly operator dependent. Here, we introduce a method for real-time quantitative stiffness mapping using an unmodified commercial endoscope, with the goal of enhancing diagnostic capabilities and restoring mechanical feedback during endoscopic procedures. Our approach combines shear wave elastography with speckle imaging and an innovative synchronization strategy that enables the measurement of shear wave propagation using an unmodified commercial endoscope. The resulting wave fields are analyzed with the noise-correlation-inspired (NCI) method[1], providing pixel-wise estimates of shear wave velocity and, consequently, quantitative maps of local tissue stiffness. The method demonstrated robust performance in both benchtop and endoscopic configurations. Validation was achieved on polymer phantoms as well as on ex vivo and in vivo biological tissues, highlighting its potential for minimally invasive biomechanical imaging and real-time tissue characterization.

Significance: Continuous, noninvasive, and depth-resolved monitoring of blood-flow-related tissue dynamics remains an important unmet need. Speckle-contrast optical spectroscopy (SCOS), including interferometric implementations such as iSCOS, provides a scalable optical route to blood-flow sensing, but conventional continuous-wave approaches lack intrinsic depth selectivity. Time-of-flight (TOF) gating offers a way to separate superficial and deeper dynamic contributions in layered tissues, such as skin-muscle or scalp-cortex, by resolving photon path lengths. Aim: We introduce a swept-source, single-channel implementation of interferometric speckle-contrast optical spectroscopy (iSCOS) to obtain TOF-resolved temporal speckle contrast, {kappa}^2, from the measured field autocorrelation g_1, and evaluate its feasibility for depth-resolved blood-flow sensing. Approach: A swept-source iNIRS system operating at 780 nm acquired interferometric signals, which were Fourier-transformed along the optical-frequency axis to recover complex TOF-resolved speckle fields. Temporal speckle contrast was then estimated at each TOF gate indirectly from g_1 using the speckle-visibility relation. Diffusion-based numerical simulations were first used to compare the direct variance-based estimator and the indirect g_1-based estimator under varying reduced scattering coefficient, diffusion coefficient, additive noise level, and bi-layer geometry. Because the simulations showed that the g_1-derived {kappa}^2 estimator was substantially less sensitive to additive noise than the direct estimator, this estimator was used for the main phantom and in vivo analyses, while the direct estimator served as a simulation comparator. The g_1-derived estimator was then applied to liquid and bi-layer phantoms, followed by proof-of-concept in vivo measurements on the human forearm during cuff occlusion and on the forehead during a Sudoku task. Results: TOF-resolved kappa2 curves recovered with the g_1-derived estimator matched DWS theory across scattering coefficients, photon path lengths, and exposure times. The estimator preserved theoretical accuracy for additive noise amplitudes up to 50% of the field amplitude, whereas the direct variance estimator showed substantial noise-induced bias and required correction. Bi-layer simulations and phantom experiments reproduced the predicted direction and onset of TOF-dependent decorrelation-rate trends in layered media. In vivo, the recovered blood-flow index tracked the expected TOF-dependent cuff-occlusion and reactive-hyperemia response in the forearm. During the single-subject Sudoku task, the left-forehead recording showed a TOF-dependent relative blood-flow-index increase of +0.8 {+/-} 1.9% at TOF = 400 ps, +9.8 {+/-} 2.2% at TOF = 600 ps, and +15.2 {+/-} 5.6% at TOF = 800 ps. This pattern is consistent with increased sensitivity to deeper tissue at longer photon path lengths, but requires cohort-level validation before quantitative interpretation as cognitive activation. Conclusions: Coupling temporal speckle-contrast analysis with swept-source iNIRS yields a proof-of-concept, depth-resolved platform for blood-flow sensing. By estimating TOF-resolved speckle contrast through the g_1-derived {kappa}^2 route, TOF-iSCOS suppresses additive-noise bias while preserving sensitivity to deeper dynamic tissue layers. The present single-channel results bridge continuous-wave iSCOS, interferometric NIRS and time-domain diffuse correlation spectroscopy (TD-DCS), and motivate future multi-channel and cohort studies for scalable cortical hemodynamic monitoring.

Three-dimensional photoacoustic tomography (3D-PAT) enables noninvasive structural and functional imaging with optical absorption contrast and ultrasonic detection depth. However, its spatial resolution is limited by acoustic diffraction, and incomplete detection geometry can substantially degrade image fidelity and quantitative accuracy. Here, we present a ULM-guided model-based reconstruction framework, termed 3D-PAULMprior that incorporates sub-diffraction vascular priors from concurrent ultrasound localization microscopy (ULM) into 3D photoacoustic reconstruction. The method uses weighted regional Laplacian regularization to integrate high-resolution vascular information into the inverse problem, thereby enhancing vascular sharpness, suppressing limited-view artifacts, and improving blood oxygen saturation estimation. We validated 3D-PAULMprior using numerical simulations, tissue-mimicking phantoms, and in vivo mouse brain imaging. Compared with conventional reconstruction, 3D- PAULMprior improved spatial resolution by over 50%, increased contrast-to-noise ratio by 261.2%, and enhanced structural similarity index by 24.6%. In vivo, 3D-PAULMprior recovered vascular structures that were poorly resolved or missing in conventional reconstructions and produced more spatially confined sO2 maps. These results establish 3D-PAULMprior as a robust multimodal reconstruction strategy for high-resolution structural and functional photoacoustic imaging.

Microbubble contrast-enhanced ultrasound (CEUS) relies on discriminating nonlinear bubble signals from linear tissue backscattering. While Singular Value Decomposition (SVD) filtering improves this discrimination, existing techniques often fail to retain the slowly-moving microbubble signals from static clutter. Here, we present a novel multi-stage singular value decomposition (MS-SVD) framework for ultrafast CEUS imaging. Our method employs plane-wave transmissions at multiple angles and acoustic pressure levels (implemented via duty-cycle modulation) and alternating transmit polarity. The beamformed data are then processed by three sequential SVD filters: (1) spatial-angular SVD to extract coherent signals across all transmit angles, (2) spatial-pressure SVD to separate linear fundamental and nonlinear harmonic components, and (3) spatiotemporal SVD to isolate moving microbubble echoes from tissue clutter. In in vitro flow phantoms and in vivo rat brain through a cranial window, MS-SVD dramatically improves microbubble detection compared to conventional SVD filtering, MS-SVD yields much stronger vascular contrast and suppresses tissue clutter to a greater extent. The resulting power-Doppler and super-resolution maps are notably cleaner and more complete: MS-SVD detects substantially more microbubble events in ULM, revealing finer vessel details and more accurate flow speeds. By capturing the full acoustic signature of microbubbles (both fundamental and harmonic), MS-SVD achieves higher contrast-to-noise and sensitivity in CEUS. These gains make it a powerful front-end for super-resolution ultrasound localization microscopy and other high-sensitivity microvascular imaging applications.

The spatial resolution of optical imaging systems is fundamentally restricted by the diffraction limit. However, in widefield live-cell microscopy, the achievable resolution is further constrained by the specimen motion, which indicates the existence of a fundamental spatio-temporal resolution trade-off between signal accumulation during the full frame integration and the resulting motion blur. To improve the fidelity with which moving objects can be imaged, a quantitative understanding of this spatio-temporal trade-off is necessary. Here, we present a systematic analysis of motion-induced resolution dynamics measured with spectral signal-to-noise ratio (SSNR). We developed a simulation framework which models the image formation of objects undergoing arbitrary motion, to evaluate the degradation of the spatial resolution under translational and rotational dynamics. Our results demonstrate that for translating objects, the spatial resolution is anisotropically reduced as a function of the orientation of the object relative to the motion vector, leading to the spectral signal-to-noise ratio degrading by up to 50% and the resolution by up to 40% for a 90{degrees} change in the motion direction. Furthermore, we show that for rotational motion, conventional radially averaged metrics such as the Fourier Ring Correlation are not able to quantify the effects of angular blur. On the other hand, the SSNR is able to accurately quantify this degradation. These findings underscore the necessity of an object-oriented imaging approach, in which acquisition parameters such as exposure time are tuned to specific biological spatio-temporal characteristics to optimize the trade-off between motion blur and spatial fidelity.

High bit-depth medical images preserve subtle intensity variations that are important for quantitative analysis and clinical interpretation, but their large dynamic range poses challenges for efficient compression. We propose a bit-plane-aware dual-stream compression framework for 16-bit medical images by separately modeling the most significant bit (MSB) and least significant bit (LSB) components. The MSB structural stream is encoded using JPEG coding with a Duplicate Segment Skipping (DSS) strategy to exploit spatial and segment-level redundancy, while the LSB detail stream is compressed using learned image compression to represent residual variations and fine-grained details. Experiments on four MRI and CT datasets show that the proposed method consistently outperforms representative traditional and learning-based codecs, achieving the lowest bit rate across all datasets. Meanwhile, it preserves high reconstruction fidelity. As a downstream application, we further demonstrate that the compressed bitstreams can be effectively integrated with DNA encoding and converted into sequences with favorable biochemical properties.

Optical microscopy is fundamental to modern life-science research, yet interpreting its results requires precise modelling of point spread functions (PSFs) within complex environments. This manuscript introduces a versatile and efficient approach to wave-optical PSF calculations that extends existing frameworks by incorporating detection PSF modelling through the principle of reciprocity. Accompanying this work is a free MATLAB software package centred on a single, minimalistic core function, PlaneWaveExc.m, which utilizes a plane-wave superposition based on the Richards-Wolf model. Despite its simplicity, the framework accounts for "real-life" complexities such as systemic aberrations, arbitrary amplitude and phase modulations, and stratified media with complex-valued refractive indices. We demonstrate the softwares broad applicability through diverse case studies, including single-molecule imaging, STED microscopy, the segmented aperture of the James Webb Space Telescope, and coherent wide-field iSCAT microscopy. Each example is supported by dedicated scripts to facilitate adaptation for specific research needs.

Light-sheet fluorescence microscopy (LSFM) has revolutionized biological imaging by enabling high spatial and temporal resolution with minimal photodamage. However, conventional LSFM techniques often suffer from striping artifacts in the resulting images due to light scattering and absorption within samples, leading to uneven illumination that negatively impacts the accuracy of subsequent image analyses. To address this limitation, we introduce dodecagon light-sheet fluorescence microscopy (dodecaLSFM), a novel approach that maximizes angular diversity to achieve homogeneous illumination and suppress striping artifacts. dodecaLSFM employs diffraction optics and cylindrical lenses to generate twelve light sheets, providing 360 degree omnidirectional illumination that significantly enhances illumination uniformity compared to traditional mSPIM, mDSLM, and ultramicroscopy systems, which use only one or two illumination planes. We demonstrate the effectiveness of dodecaLSFM by achieving high-resolution imaging of whole mouse brain vasculature following tissue clearing, allowing precise morphometric analysis of vascular networks without striping artifacts. Furthermore, we show that combining dodecaLSFM with expansion microscopy (ExM) enables whole-organ 3D imaging at cellular resolution. This novel approach provides an advanced, scalable solution for large-volume imaging, facilitating detailed structural and functional studies across diverse biological applications.

Three-dimensional (3D) imaging represents the development of next generation of fluorescence microscopy. However, routine axial down-sampling makes isotropic resolution unrealistic. Here, we propose DeepUI, a physical zero-shot framework designed to achieve isotropic 3D fluorescence images from a low axial sampling rate. DeepUI fully leverages the intrinsic characteristics of 3D images through physics-guided degradation, which incorporates spatial-frequency joint learning to generate a scaled optical transfer function, combined with noise degradation and an up-sampling branch. Typically requiring just 5 minutes for training and 0.5 minutes for high-throughput and fast prediction, we demonstrate the superior performance of DeepUI to get isotropic results, and the exclusivity to axial down-sampling conditions, even in more challenging conditions, including defocused background, noise, and resolution blur.

Optical fluorescence microscopy enables visualization of biological structures and dynamics. However, the intrinsic diffraction limit, especially axially, and depth-related scattering noise compromise the image resolution and fidelity. Computational 3D deconvolution is a promising approach for mitigating these issues, yet its execution is hindered by inaccurate and cumbersome theoretical modeling or experimental measurement of 3D point spread function (PSF), as well as ineffective 3D noise regularization. Furthermore, in the 3D super-resolution regime, there remains a lack of standardized tools for evaluating 3D super-resolution fidelity. Here, we present the 3D adaptive deconvolution and evaluation (3D-ADE) toolkit, which comprises 3D-Ada deconvolution with physics-oriented automatic 3D-PSF calibration, and 3D-SQUIRREL for 3D super-resolution quality assessment. It effectively resolves noise instability, eliminates the need for 3D-PSF calibration, and reliably assesses the fidelity of 3D resolution extension via deconvolution, physical, and deep-learning-based methods. Accessible via multiple software platforms, 3D-ADE enhances the versatility of 3D deconvolution and fills the gap in 3D super-resolution evaluation tools, and thereby advances volumetric fluorescence imaging applications.

We introduce a workflow to identify oligomeric structures that are recorded with single-molecule localization microscopy (SMLM) under cryogenic conditions. Typically, these oligomers are assumed to consist of protomers arranged as equilateral two-dimensional polygons and every protomer is labeled with a dye molecule for visualization. Unlike previous work, we consider scenarios in which the sample plane has an unknown orientation relative to the focal plane. Our contribution is a high-precision plane-fitting algorithm to determine the sample plane, combined with geometrical transformations and two circle-fitting algorithms to identify the oligomeric structures. Our simulations on synthetic data demonstrate that the proposed workflow achieves high accuracy in estimating both the unknown tilted plane and the oligomer size.

Point-of-care ultrasound has transformed bedside diagnostics, yet current systems remain limited by rigid form factors, bulky external electronics and the need for skilled operators. Here we report a conformable ultrasound imaging patch that integrates a 1024-channel CMOS ultrasound application-specific integrated circuit (ASIC) directly beneath a conformable piezocomposite transducer array. The 10 mm X 8 mm, 1024-element ASIC contains on-chip transmit and receive beamforming, reducing the effective off-chip channel count by 16X while preserving image fidelity. Fabricated on a flexible polyimide substrate and bonded using anisotropic conductive film, the patch operates untethered from conventional ultrasound consoles and requires only a laptop for control and data acquisition. The device supports focused, plane-wave and diverging-wave transmission with steering over {+/-}30{degrees} in azimuth and {+/-}15{degrees} in elevation, achieving peak-to-peak acoustic pressures up to 7 MPa at a 4.4-MHz center frequency (mechanical index of 1.7), within diagnostic safety limits. Phantom experiments demonstrate three-dimensional imaging with axial and lateral resolutions (in both XZ and YZ planes) of 0.5 mm and 2 mm, respectively, and accurate contrast reproduction in tissue-mimicking phantoms. Human studies further demonstrate three-dimensional (3D) visualization of the internal jugular vein and carotid artery, as well as rib-shadow-free imaging of pleural motion during respiration. This work establishes a scalable architecture for chronic, wearable ultrasound imaging and highlights the potential of CMOS-integrated, conformable ultrasound systems for continuous physiological monitoring and remote diagnostics.

Coherent structured illumination microscopy (c-SIM) is a synthetic aperture optical technique for sub-diffraction limit imaging that extends the utility of traditional SIM to non-fluorescent samples. Here, we present a complementary 5-beam implementation of c-SIM that provides enhanced optical sectioning compared to conventional quadrupolar illumination. Since our approach detects intensity images due to coherent light scattering, it avoids the complications associated with detecting complex fields. Through comparative measurements on calibration samples and live microalgae, we show that 5-beam c-SIM effectively suppresses coherent defocus effects, improving image quality while simultaneously providing a 2-fold lateral resolution improvement.

Objective: This study aims to develop and validate a quantitative, depth-resolved anisotropy imaging framework that extends ARFI-based focal degree-of-anisotropy (DoA) estimation into two-dimensional mapping by modeling the depth-dependent relationship between shear modulus ratio (SMR) and peak displacement ratio (PDR). Methods: We propose APRIL (Adaptive Polynomial Regression for anisotropy Imaging via ARFI-induced DispLacements), a framework for quantitative, depth-resolved DoA imaging that adaptively selects polynomial regression or shape-preserving spline interpolation based on excitation PSF asymmetry. Training data were generated using an LS-DYNA3D + Field II simulation pipeline in homogeneous transversely isotropic media (SMR 0.9-4.9). Testing included shifted SMRs under varied acoustic conditions and three heterogeneous inclusion configurations (anisotropic inclusion in isotropic background and vice versa). Experimental validation was performed in an in-vivo murine tumor model over the time, ex-vivo chicken breast, and tissue-mimicking gelatin phantoms, using a Verasonics system with an L11-5v transducer. Results: APRIL achieved depth-resolved SMR prediction errors below 9% over 10-30 mm, with highest accuracy in the focal region (MAE 2.3%, RMSE < 0.1) and stable performance across PSF transition zones. In heterogeneous phantoms, it reconstructed anisotropy maps with SSIM up to 86% and MPE below 7%, accurately delineating inclusion boundaries. Under acoustic parameter variations, mean absolute errors remained below 10%, demonstrating robustness to system and tissue heterogeneity. Conclusion: APRIL enables robust, two-dimensional anisotropy imaging beyond focal estimates. Significance: The method provides a physically grounded and generalizable framework for clinically viable anisotropy biomarkers in muscle, tendon, kidney, tumor and breast tissues.

Three-dimensional fluorescence microscopy often exhibits anisotropic resolution because axial information is poorly sampled and more blurred than lateral information, which complicates quantitative interpretation of fine 3D structures. Although optical remedies and computational restoration have been explored, many approaches require demanding system calibration or rely on accurate PSF models and assumptions that are difficult to satisfy across all samples and modalities. Here we present DeepIso, a self-supervised isotropy restoration framework that couples supervised pretraining with an internal-learning inference stage to estimate degradation directly from the measured volume. Without explicit PSF specification or enforced lateral-axial structural equivalence, DeepIso recovers axial frequency content and improves the continuity of elongated structures while retaining fine features, with superior performance over existing computational approaches in terms of both visual inspection and quantitative metrics. The method is validated on synthetic benchmarks and experimental datasets, demonstrating isotropy enhancement across confocal, light-sheet, and 3D structured illumination microscopy, thereby supporting downstream volumetric analysis including segmentation and tracking.

For decades, the photon counting histogram (PCH) was used as the sole method to quantify fluorophore numbers in a diffraction-limited focal volume. This technique combines spatial excitation profiles, and the distribution of photon counts to register the photon emission statistics of individual fluorophores. However, this approach has not yet been transferred to widefield fluorescent imaging due to the lack of fast and single photon sensitive camera sensors which can capture the photon emission statistics of a single fluorophore. Here, we explore avenues towards quantitative analysis of the active fluorophore number by leveraging recent advancements in single photon avalanche diode (SPAD) array technology. Binary exposures of a SPAD array can be synchronized with picosecond laser pulses to measure the PCH in a widefield setting. Then, by modeling the statistical relationship between the active fluorophore number and the PCH in a region of interest following a laser pulse, we can perform Bayesian inference of this number. The model is demonstrated experimentally by counting quantum dots and various numbers of fluorescent dye molecules bound to DNA origamis. We find that this method has several important applications in widefield microscopy, including enhanced localization microscopy and constrained fitting of multiple unresolvable fluorescent emitters.